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Information technology has revolutionized the traditional structure of markets. The removal of geographical 
and time constraints has fostered the growth of online auction markets, which now include millions of economic 
agents worldwide and annual transaction volumes in the billions of dollars. Here, we analyze bid histories 
of a little studied type of online auctions — lowest unique bid auctions. Similarly to what has been reported 
for foraging animals searching for scarce food, we find that agents adopt Levy flight search strategies in their 
exploration of "bid space". The Levy regime, which is characterized by a power-law decaying probability 
distribution of step lengths, holds over nearly three orders of magnitude. We develop a quantitative model for 
lowest unique bid online auctions that reveals that agents use nearly optimal bidding strategies. However, agents 
participating in these auctions do not optimize their financial gain. Indeed, as long as there are many auction 
participants, a rational profit optimizing agent would choose not to participate in these auction markets. 



Introduction 

Animals searching for scarce food resources display move- 
ment patterns that can be statistically classified as Levy 
flights lHJ^Sl. Levy flights | 9 | represent the best strategy that 
can be adopted by a searcher looking for a scarce resource 
in an unknown environment 1 10|, and foraging animals seem 
therefore to have learned the best strategy for survival. Levy 
flights describe also the movement patterns of humans in real 
space (TlJ and the variability of economic indices [12J, but 
these observations do not correspond to search processes as 
in the case of foraging animals. Surprisingly, there is no in- 
dication of whether humans also use Levy flight strategies 
when searching for scarce resources. Analyzing apparently 
unrelated data regarding online auctions, we address here this 
question and show that, when searching for scarce resources, 
humans explore the relevant space in the same class of strate- 
gies as foraging animals do. 

Lowest unique bid auctions are a new generation of online 
markets 1 13 -181 . Agents winning lowest unique bid auctions 
may purchase expensive goods for absurdly low prices; cars, 
boats and even houses can be bought for only hundreds of 
dollars. The idea of the auction is strikingly simple. A good, 
typically with a market value V of at least a thousand dollars, 
is put up for auction. The auction duration is fixed a priori. 
A bid can be any amount from one cent to a pre-determined 
maximum value M, generally lower than one hundred dollars. 
Each time an agent makes a bid on a value 1 < b < M, she 
pays a fee c, which ranges from one to ten dollars depending 
on the auction. During the bidding period, an agent knows 
only the status of her new bid, that is, whether it is winning 
or not. None of the agents knows on what values the other 
agents have bid until the end of the auction. When the bidding 
period expires, the agent who made the lowest unmatched bid 
can purchase the good for the value of the winning bid (see 
Fig. 1 for an illustration of the determination of the winning 



bid). 
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Figure 1 : Unique bid auctions. Illustration of the rules of a lowest 
unique bid auction. At the end of the auction, the winner results 
to be the agent who has bid 30, which represents the lowest unique 
bid. All other bids are not unique apart from the one of 50, which is 
not the lowest one. In highest unique bid auctions the mechanism is 
reversed, and the winner is the agent making the highest unmatched 
bid. Illustration of the rules of a lowest unique bid auction. At the 
end of the auction, the winner results to be the agent who has bid 30, 
which represents the lowest unique bid. All other bids are not unique 
apart from the one of 50, which is not the lowest one. In highest 
unique bid auctions the mechanism is reversed, and the winner is the 
agent making the highest unmatched bid. 



Lowest unique bid auction markets are competitive arenas. 
Each agent performs a search for a single target whose po- 
sition changes from auction to auction, as it is determined by 
the bid history of the whole population of agents. Since the 
cost of each bid is as much as 100 times larger than the natu- 
ral unit of the bid, the number of bids that can be made by a 
single agent is limited and allows only a partial exploration of 
the bid space. Successful agents need to identify good strate- 
gies in order to maximize their winning chances and thus limit 
their risk. 
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Figure 2: Individual activity. (A) Bid values explored by agent 
1632 on auction 19 in the data set www . uniquebidhomes . com. 
Bids are sorted chronologically, and the figure reports the value bt 
of the t-th bid. The unit of the bid amount is one hundredth of 
an Australian dollar. (B) Absolute value of the difference between 
two consecutive bids. The exploration of the bid space is character- 
ized by a bursty behavior, where many small movements are occa- 
sionally followed by large jumps. (C) Cumulative distribution func- 
tion of the change in bid value. The distribution is well fitted by 
a power-law, with decay exponent consistent with a = 1.3 ± 0.1 
(dashed line). The agent therefore explores the bid space using a 
Levy flight strategy. Notice that the curve bends down because of 
the finiteness of the bid space. (D) Probability density function 
of the Levy-flight exponents adopted by agents in lowest unique 
bid auctions (www.uniquebidhomes.com). The blue line in- 
dicates the average value (a) ^ 1.26 of the distribution, the red 
line identifies the mode at ^ 1.21 of the distribution, the orange 
lines bound the region within one standard deviation a 0.23 
from the average. (E) Probability density function of the Levy- 
flight exponents adopted by agents in highest unique bid auctions 
(www.bidmadness.com.au). In this case we find (a) ^ 1.36, 
at ^ 1.35 and a :^ 0.23. 



Lowest unique bid auctions are just a particular variant of on- 
line pay-to-bid auctions, but other types of pay-to-bid auc- 
tions are regularly hosted on the web. For example, in high- 
est unique bid auction the mechanism of lowest unique bid 
auction is inverted, and the winning bid is determined by the 
highest value closest to a pre-determined upper bound value. 
Since these auctions still involve a blind search of the win- 
ning value, highest unique bid auctions are equivalent to low- 
est unique bid auctions. Indeed, in this paper we analyze data 
taken from both types of auctions. 

Other online pay-to-bid auctions, however, can be very dif- 
ferent from lowest unique bid auctions. For example, the so- 
called penny auctions, which have acquired a great popularity 
in recent years, appear quite similar to but are not. As in the 
case of lowest unique bid auctions, the cost of the fee is at least 
100 times larger of the bid increment, and as a consequence. 



the final value of the winning bid is much lower than the real 
value of the good up for auction. However, in penny auc- 
tions the value of the winning bid is publicly known and can 
only grow during the auction (i.e., the word "penny" is used 
because, in penny auctions, bid increments are equal to one 
cent). While escalation plays a very important role in penny 
auctions, in this type of auctions agents do not need to explore 
the bid space because the value of the winning bid is known. 
Penny auctions have been the focus of some theoretical and 
empirical studies |[T9ti23]| . 



Results 

We collected data from three distinct web sites hosting low- 
est unique bid auctions. We automatically downloaded and 
parsed the content of the tables reporting the bid history of 
closed auctions. These data sets contain all the information 
on individual auctions, including the details of each bid: its 
value, when it was made and who placed it. These data allow 
us to keep track of all the movements performed on bid space 
by a given agent bidding in a specific auction. 
We show in Figure 2A a typical exploration of the bid space 
performed by a single agent. The exploration of the bid space 
is bursty: consecutive bid values are generally close to each 
other, but from time to time the agent performs "long jumps" 
in bid space. We first compute the jump lengths (Fig. 2B) and 
estimate their probability distribution function (Fig. 2C). We 
find a strikingly robust power-law scaling consistent with the 
exploration of the bid space using a Levy flight search strat- 
egy |9|. Note that here we use the notion of discrete Levy 
flights. Time and space are in fact discrete, and the exploration 
of the bid space is modeled as a discrete time Markov chain 
[with transition probability defined in Eq. [8| . Our discrete 
model converges to a standard Levy flight only in the contin- 
uum limit of space and time \24\. The power-law scaling can 
be observed both at the level of single agents (whenever the 
number of bids is sufficiently large for estimating the distribu- 
tion; c.f. Figs. 2C and Supporting Information) and globally, 
by aggregating the length of the jumps made by all agents in 
all auctions (Figs. 3 A and Supporting Information). The den- 
sity distribution of the exponents calculated over single agents 
is peaked around a mean value (a) 1.3 (Figs. 2D, 2E and 
Supporting Information), the same exponent value we esti- 
mate for the aggregated data. Significant variations around 
the average value are anyway present, and reflect the hetero- 
geneity of the agent strategies. The density distributions of 
Figs. 2D and 2E are in fact calculated by considering different 
agents bidding in different auctions. 

The power-law scaling and its measured exponent are very 
stable. Exponent estimates do not depend on the direction of 
the jumps (Figs. 3B and Supporting Information) or the level 
of activity of the agent (Figs. 3C and Supporting Information). 
Surprisingly, performing Levy flights does not appear to be 
a learned strategy. Instead it appears to be an intrinsic fea- 
ture of the mental search process: the jump lengths in the bid 
space follow the same power-law at any stage of the auction 
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(Figs. 3D and Supporting Information). 

Our results represent the strongest empirical evidence for the 
use of Levy flight strategies in the search of scarce resources 
reported in literature up to now. Differently from previous 
studies where "two orders of magnitude of scaling can repre- 
sent a luxury" |6|, here the power-law decay can be clearly 
observed even over four orders of magnitude. It is unlikely, 
though, that adopting Levy flight strategies is a deliberate 
choice of the agents, just as it is not likely that animals search- 
ing for food consciously follow a Levy flight strategy. Nev- 
ertheless, the data demonstrate that the changes in bid value 
are statistically consistent with a power-law decaying distribu- 
tion over several orders of magnitude (see and Supporting In- 
formation) |25|. Simple correlation measurements show also 
that the lengths of consecutive jumps are independent of each 
other (see and Supporting Information). We believe that the 
power-law is valid over such a broad regime because the space 
is not strictly physical. That is, movements of tens of thou- 
sands of cents can be performed for the same cost of those of 
only one cent. Agents thus explore the bid space in an effec- 
tively super-diffusive fashion, and steps are made with infinite 
velocity. 



Model 

Next, we model the lowest unique bid auction process. Con- 
sider N agents competing in a lowest unique bid auction. 
We model the successive bids of these agents as Levy flight 
searches on bid space. Each agent moves in a bounded one- 
dimensional lattice with an a priori chosen exponent value, 
which may be regarded as the agent's strategy in the auction. 
In our formulation, every agent performs the same number T 
of bids and may return to already visited sites. At the begin- 
ning of the auction, every agent sits at the leftmost site on the 
lattice and then performs T movements by changing, at each 
step, her actual position by an amount randomly drawn from a 
power-law distribution. If at stage t — 1 the agent with strategy 
a is sitting at position j , then at stage t she jumps to position 
i with probability proportional to |i — j |~^. This model pro- 
vides us with an independent way to determine the exponent 
values of the Levy flights and offers a strikingly good statisti- 
cal description of the data (Fig. 2B and Supporting Informa- 
tion). 

We focus our attention on a generic agent bidding with strat- 
egy 13 and on her chances to win auctions in which the rest of 
the population is bidding with strategy 7. More complicated 
situations may in principle be studied with the same formal- 
ism. 
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Figure 3: Bidding strategies of agents are Levy flights. (A) 

Probability density function of the bid change for all agents 
in all auctions. We analyze data sets from three differ- 
ent web sites hosting auctions: www.uniquebidhomes.com 
(black circles), www.lowbids.com.au (red squares ) and 
www.bidmadness.com.au (blue diamonds). (B) Probability 
density function of positive (black circles) and negative (red squares) 
bid changes. (C) Probability density function of the change amount 
for data aggregated over agents with different levels of activity (T in- 
dicates the total number of bids made by an agent in a single auction). 
(D) Probability density function of the change amount at different 
stages of the auctions (t stands for order of the bid change in the bid 
history of an agent). In (A), (B) and (D) results have been obtained 
for lowest unique bid auctions (www.uniquebidhomes.com). 
All dashed lines stand for best power-law fits (least square) and all 
exponent values are consistent with a = 1.4 ± 0.1. The unit of the 
bid value change amount is one hundredth of an Australian dollar. 



bids on value i is 



7m (7) 



(1) 



with m (7) = Y^^i proper normalization constant. Here 
we consider the simple case in which all agents adopt the same 
bidding strategy 7. The probability of Eq. ([T]) can be anyway 
made more general by assuming that agents chose strategies 
from a density distribution g {a) and calculating the proba- 
bility of Eq. lias p (i) = J da i~^/m {a) g {a). After all 
agents have bid, there will be rik bids on the k-th bid value. 
Such variables clearly obey the constraint N = Xl^i rik. 
The probability to observe a particular configuration {n} = 
(ni, 712, . . . , n/e, . . . , tim) is given by 



Single bid 

Consider first the case in which agents make a single bid. The 
probability that a generic opponent, using bidding strategy 7, 



M 



P, ({n}) = iV! n 



/C = l 



(2) 



which is a multinomial distribution with weights given by 
Eq. ([T]). In particular, the probability that only one bid (i.e.. 
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Figure 4: Model predictions. Economic return r/3,7 [Eq. (|7])], di- 
vided by the number of bids T, of an agent bidding with strategy /3 
when competing, in a lowest unique bid auction with upper-bound 
M = 1, 000 and for a good of value V = 10, 000, against N = 100 
opponents bidding with strategy 7. Unless specified, the quantity 
reported in this plots is computed by numerically solving the 
equations of the model. (A) Case where each agent performs a sin- 
gle bid in the auction, for three values of 7. Theoretical predictions 
(lines) are compared with the results of numerical simulations (sym- 
bols). In each simulation of the auction, we randomly extracted N 
bid values j with probability proportional to , and a single bid 
value V with probability proportional to . For a given set of 
parameters, we repeated the same simulation G — 10,000 times, 
and calculate the number of times g in which the bid value extracted 
from the power-law distribution with exponent (3 was the winning 
bid, and the sum / of these winning bid values. The economic re- 
turn has been finally calculated as r/3,7 = {gV — I) /G. (B) Explo- 
ration of parameter space reveals the existence of a saddle point at 
/3s = 7s ^ 1.27. (C) Case where each agent performs 10 bids in the 
auction, for three values of 7. Numerical simulations have been car- 
ried out as in the former case, but considering agents moving in the 
bid space according to Eq. |8]). (D) Exploration of parameter space 
reveals the existence of a saddle point at /3s = 7s — 1-38. 



a unique bid) is made on value i is 



^7 (0 Pi i^i = 1) = E 



Np^ (i) {i)] 



N-l 



-1^7(W) 



(3) 



Focus now on the agent with bidding strategy p. The proba- 
bility that, making a bid on value v, she makes a lowest unique 
bid can be calculated exactly by summing the multinomial dis- 
tribution of Eq. ([2]) over all configurations for which there are 
no bids on the value v and there is not a unique bid on a value 
smaller than and finally multiplying this factor by the prob- 
ability that the agent with bidding strategy {3 bids on the value 
V. Such exact calculation is however unfeasible due to the ex- 
tremely high number of possible combinations, and therefore 



we approximate the probability that, making a bid on value v, 
the agent with bidding strategy j3 makes a lowest unique bid 
as 

l^^^ {v) = {v) [1 - p^ {v)f \[[l-u^ m . (4) 

k<v 

The r.h.s. of Eq. ^ is the product of three terms: p/3 (v) is 
the probability that the agent bids on value [1 — p^ {v)]^ 
is the probability that none of the opponents have bid on 
value v; Yli^^y [1 — {k)] is the probability that none of 
the bid values smaller than v are occupied by a single bid 
made by one of the opponents. In spite of the fact that 
Eq. ^ is just an approximation of the real /^^^ {v), the ap- 
proximation can be considered good because able to repro- 
duce the results obtained from the direct simulation of the 
process (see the section Results). Moreover in the simplest 
case in which TV = 1, it correctly reduces to the exact value 

hn (^) = (^) Uk<v [1 - Pi l^)]' 

Finally, the probability that the agent with bidding strategy (3 
wins the auction is 



M 



^/5,7 = Yl ^/5,7 (^) 



and, on average, the value of her winning bid is 



M 



(5) 



(6) 



Repeated auctions 

Imagine now to repeat the same auction G independent times. 
The probability that the agent bidding with strategy (3 wins g 
times out of G total auctions is given by a binomial distribu- 
tion 



^/5,7 (9) 



G-g 



If the agent with bidding strategy (3 wins g auctions, the sum 
of her winning bids is a random variable / whose probability 
is determined by 

Vi-\-V2-\-...-\-Vg=I 

where the sum runs over the integer indices Vi, V2, . . . , Vg 
with the constraint that their sum should equal /. Excluding 
bidding costs, the average return of the agent in g victories is 

rp,^(9) = (9V -I)/G . 

In general, the probability that the sum of the winning bids 
is equal to / in an arbitrary number of auctions won by the 
player with bidding strategy j3 can be calculated as 

Rp,^{I) = Y^Pp,^{g) Rp,^{I\g) , 
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and a similar expression can be derived for the distribution 
of {g). However, we are interested in the case in which 
the number of auctions diverges (G 1). In this Hmit, 
we can approximate the number of victories with its aver- 
{9) — Gwj3^^ as well as the sum of the winning bids 
as / = {g){v)i3^j = Gw^s^j (^)/3,7- The return of the agent 
with bidding strategy (3 is therefore 



(7) 



For > c, the agent has a positive return for participating 
in the auction, whereas, for < c, her return is negative. 



Multiple bids 



Given a generic agent with bidding strategy a, her first bid is 



Ym (a). For 



placed on value i with probability qa'' {i) 
the subsequent bids, we need to define a transition matrix Qa, 
whose generic element {Qa)ji gives the probability that the 
agent bids on value i when her previous bid has been made on 
value j. In our model, we have 



rrij {a) 



(8) 



for all i and j in the interval [1, M]. 5 {•) is the Kronecker 
delta, equal to one if its argument is equal to zero, and equal 
to zero otherwise. The normalization constant ruj {a) = 

YlfLi i^j ~ ensures the proper definition of the tran- 
sition matrix. The matrix Q describes a random walker per- 
forming uncorrected Levy flights with exponent a. Notice 
that the agent has no memory of her previous bid values and 
therefore she may place more than a bid on the same value. At 
the generic step t, the probability that the agent with bidding 
strategy a bids on the value i is 

M 

The probability that this agent has bid, during her T bids, on 
value i is then 



The term 1 — qa"^ (i) counts the probability that the agent has 
not bid on value i at stage t. The probability that the agent has 
not bid on value i at any stage is therefore the product of this 
single step probabilities. Finally, the probability that the agent 
has bid on value i at least once is calculated as the probability 
to have bid on value i an arbitrary number of times minus the 
probability to have never bid on value i. 

Now go back to the situation in which an agent with bidding 
strategy (3 is opposed to a population of N agents with bidding 
strategy 7. The probability that the agent with bidding strat- 
egy /3 has bid, in T steps, at least once on value i is 5^^^ (z). 



The probability that one of the opponents, bidding with 
strategy 7, makes a unique bid on value i is given by 



4^(i)=iVsf) (i) 



1 



lN-1 r 



AT) 



,iT) 



(0 
(T) 



the 



1 



ST) 



(0 



product 

AT-l 



of 



two 



(9) 
terms: 



is the probability that 



bid on value i is unmatched by any of the other — 1 
opponents, while 1 — s^^^ (i) is the probability that also the 
agent, with bidding strategy (3, does not bid on value i. The 
probability that the agent with strategy /3 wins the auction 
with a bid on value v is 



ST) 



{v) 



1 



ST) 



(v) 



N 



XT) 
^/3,7 



(V) 



k<v 



(10) 

respectively standing for the product of the probabilities that: 
she bids on value v, none of the other agents bids on value v; 
none of the bids with value smaller than v is unique. Eqs. ^ 
and ([To]) represent the generalization of Eqs. Q and (|4]), re- 
spectively. In Eq. ([To]) we made the same type of approxima- 
tion as the one used for writing Eq. Ui. The probability w^^}^ 
that the agent with bidding strategy p wins the auction and 
the average value (^)^^^ of her winning bids can be respec- 
tively calculated using Eqs. ([5]) and ([6]). Finally, excluding 
bidding costs, the return of the agent with strategy (3 over 
an infinite number of auctions is again given by Eq. ([T]). For 
> T c, the agent has a positive return for participating in 
the auction, whereas, for < T c, her return is negative. 



Model predictions 

We show in Fig. 4 the results obtained with our analytical 
model. The presence of a saddle point at 7^ = indicates 
that is an optimal strategy or Nash equilibrium 126^281 . 
When the opponents do not bid rationally (i.e., 7 7^ 7^), it is 
more convenient to use a strategy /3 7^ /3s . On the other hand, 
when the other agents bid rationally (i.e., 7 = 7s), there is no 
better strategy than Ps- The value of depends on the pa- 
rameters A" and T , but for realistic choices (see and Support- 
ing Information and Fig. 4), /3s is in the range 1.2 to 1.5, the 
same range of the exponent values we estimated from the data. 
Thus, despite its simplicity, our model captures the main fea- 
tures of the real auctions. Performing Levy flights with small 
exponents (ballistic motion) yields unique bids that are un- 
likely to be the lowest. On the other hand, performing prefer- 
entially short jumps (high exponents, diffusive motion) guar- 
antees to always bid on small values which are unlikely to be 
unique. Intermediate values of the exponent (super-diffusive 
motion) represent a compromise between staying low and be- 
ing unique, and therefore lead to maximal winning chances. 
These considerations are valid only for finite values of A^ and 
T, which is the realistic case. Because the available positions 
in the lattice are finite, when either N orT grow, the probabil- 
ity to observe a unique bid progressively approaches zero L29 J . 
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Figure 5: Economic return of agents. Auction winners 
tend to pay half of the value of the good, while auctioneers 
tend to earn twice the value of the good. We consider two 
different data sets, one regarding lowest unique bid auctions 
[www . uniquebidhomes . com, (A) and (C)] and the other highest 
unique bid auctions [www . bidmadness . com. au, (B) and (D)]. 
The unit is one hundredth of an Australian dollar. (A) Relation be- 
tween the relative cost (Tu, c + 6*) to the winner of an auction 
and her income V. Tu, indicates the number of bids made by the 
winner of the auction, c is the cost of the fee, b* indicates the value 
of the winning bid, and V is the value of the good put up for auc- 
tion. Each point represents an auction. The gray area corresponds to 
the region of negative return for winners of the auctions. The black 
line indicates the average value of the relative cost to the winners. 
Auctions are, on average, very profitable for the agent winning the 
auction, but the probability of a specific agent winning the auction is 
very low. Data refer to the data set www . uniquebidhomes . com. 
(B) Relative income (B c + 6*) of the auctioneers as a function 
of their investment V. B is the number of bids placed by all agents. 
The gray area denotes the region of negative return for the auction- 
eers. The black line indicates the average value of the relative profit 
of the auctioneers. On average, organizing auctions is very profitable. 
Data refer to the data set www . uniquebidhomes . com. (C) Re- 
lation between the relative cost to the winner of an auction and her 
income for the data set www . bidmadnes s . com . au. (D) Relative 
income of the auctioneers as a function of their investment for the 
data set www . bidmadness . com . au. 



Notice that at the saddle point 7^ = f3s, slU N -\- 1 agents are 
using the same bidding strategy and therefore they all have the 
same chances to win the auction. In particular, the probability 
that a generic agent wins the auction is 1^/3^,7^ < 1/ (A^ + 1), 
where the inequality may arise because a unique and lowest 
bid may not exist. 

The value of the exponent, corresponding to the optimal Levy 
flight strategy in lowest unique bid auctions, is distinct from 
the one found in the case of purely random searches 1 10], and 
empirically observed in the movement patterns of foraging an- 



imals |T-^. The quantitative difference arises, we believe, as 
a consequence of the anisotropy of the bid space (low values 
are favored), the role of competition, and, more importantly, 
the fact that the target is not "static" but moving according to 
the actions of the whole population of agents. 



Discussion 

In lowest unique bid auctions, agents have the possibility to 
win goods of high value for impossibly low prices (Figs. 5 A 
and 5C), However, these all-pay auction markets are designed 
to be very profitable for the auctioneers 1 30 - 33 1, who, on aver- 
age, double their investment (Figs. 5B, 5D and Supporting In- 
formation). For auctioneers, the profitability of lowest unique 
bid auctions is in fact guaranteed by the validity of the in- 
equality V < B c, where B stands for the total number of 
bids and equals (TV + 1) T in our model. Under this con- 
straint however, the payoff of a generic agent in a perfectly 
rational population is always negative since 



V 



7V + 1 



<Tc 



and there is no expected economic gain to be obtained for par- 
ticipating as a bidder in the auction markets. The rationality 
of the economic agents in adopting optimal strategies seems, 
therefore, in contrast with the ultimate irrationality that in- 
duces agents to take part in these auction markets. 

Competitive irrationality, based on rational choices, has been 
investigated in economic theories |34-37|, such as the dollar 
auction game (381 . The decision to participate or not partici- 
pate in lowest unique bid auctions presents a paradox for po- 
tential bidders. If the number of agents participating in the 
auction is not too high, then the auction would bring a posi- 
tive economic return to the agents, but not to the auctioneers. 
For example, in the case in which only one bidder participates 
in the auction, this bidder would have the maximal economic 
return by placing a single bid on the lowest value allowed. But 
by this token, every agent will feel that participating is prof- 
itable as long as not many other agents have bid yet. However, 
no agent can know how many other agents will actually bid on 
the good. 

Our results raise a number of important research questions. 
First, which brain regions are responsible for implementing 
the search strategies used by agents? Since agents use similar 
search strategies to bees or birds, it is likely that there is no 
frontal cortex involvement. Using neuroimaging techniques 
such as fMRI it should be possible to answer this question. 
Second, does the economic paradox that the agents face reveal 
itself in brain activity patterns? Specifically, do some of the 
changes in brain activity observed for preference reversal 
l40l occur also in this case? Additionally, our results suggest 
that controlled lowest unique bid auction markets would offer 
the possibility to run large-scale experiments at relatively low 
cost II4TII . These experiments could be used for monitoring 
the behavior of agents in auction markets with tunable optimal 



search strategies, and see if (and how fast) agents are able to a version of these data at the 

adapt their behavior to optimaHty. web page f ilrad . homelinux . org/resources. 
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S 1. "LOWEST UNIQUE BID" AND "HIGHEST UNIQUE BID" AUCTIONS 
A. Description of the auction 




Bx c/V 

Figure SI: Cumulative distribution of the relative return (B x c) /V for the auctioneers. For www . uniquebidhomes . com data set (black 
line) the 95% of the auctions are in the gray region of positive return, where the relative return {B x c) /V > 1. On average, the relative 
return of the auctioneers is 2.2. For www . bidmadne s s . com . au data set (red line) the 72% of the auctions have produced a positive relative 
return and on average the relative return is 1.6. 



Lowest Unique Bid (LUB) auctions are special on-line auctions which have reached a considerable success during last years. 
Their peculiarity consists in the fact that they are reverse auctions: rather than the bidder with the highest bid (as in the case of 
traditional auctions), the winner is the person who makes the LUB (see Figure la in the main text). The rules of a LUB auction 
are very simple. At the beginning of each auction, the auctioneers put up for auction a good of value V. After the beginning of 
the auction and for a certain period of time (in general of order of weeks), agents participate to the auction by making bids. The 
natural unit of the auction is one hundredth of dollars, euros, etc (i.e., the currency depends on the country where the auction is 
hosted). Bids may be any amount (in cents) between one cent and a maximal bid amount M (generally lower than ten hundreds 
of cents). Sometimes, the value of M is not fixed, but the bid space is anyway naturally bounded since none of the agents wants 
to bid more than V. The value of V depends on the auction, but generally its order of magnitude is of thousands of hundreds 
of cents. Making a bid costs a fee c (typically from one hundred to ten hundreds cents). After each bid b, the agent receives an 
automatic message, from the web site hosting the auction, saying whether that bid was the winning bid (i.e., the bid is the LUB) 
or not. The agent is constantly informed about the status of her bids (i.e., whether one of them becomes the LUB or is not longer 
the LUB). However, each agent knows only what she has bid, without any information on which values the other agents have 
bid. In general, there is no restriction for the number of bids that the same agent may place. When the time dedicated to the 
auction expires, the winner is the agent who made the LUB and can therefore purchase the good for the value of her winning bid. 
If at the close of the auction a single lowest unique bid does not exist, the successful bid becomes the lowest one made by only 
two agents and the winner is the one who has bid first on such value. In the case in which also a bid made by only two agents 
does not exists, then the winning bid becomes the lowest one made by only three agents and so on. Again in these situations, 
the winner is the agent who has first placed a bid on the winning value. In our data sets however, we always observe that the 
winning bid is an unmatched bid. 

There are several slight variations of this kind of auctions. Very often, the end of the auction is not determined by an expiration 
time, but by a minimum required number of bids, a priori fixed by the auctioneers. In other variations called Highest Unique 
Bid (HUB) auctions, the winning bid is the unique one closest to M. 

Independently on the type of auctions, this kind of auctions are particularly profitable for both the auctioneers and the winners 
of the auctions. Figure 5 of the main text and Figure [ST] clearly show that there are only few exceptions in which auctioneers or 
winners have lost money, but in the majority of the auctions their returns are positive. 
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Data set 


Tot. Auctions 


Tot. Agents 


Tot. Bids 


(M) 


(c) 


(N) 




UBH 


189 


3 740 


55 041 


362 


437 


50 


6 


LB 


55 


445 


3 740 


1284 


478 


13 


6 


BM 


336 


3 719 


127 275 


504 


174 


40 


14 



Table S 1 : Summary table of the data sets analyzed in this paper. We report, from left to right, the name of the data set, the total number of 
auctions, the total number of different agents, the total number of bids, the average value of the maximal bid value, the average amount of the 
fee, the average number of agents involved in an auction and the average number of bids made by a single agent in a single auction. The unit 
of the bid values is one hundredth of an Australian dollar. 



B. Description of the data sets 



We collected data from the web sites www . uniquebidhomes . com (UBH), 
www . lowbids . com . au (LB) both hosting LUB auctions and 

from www .bidmadness . com. au (BM) organizing HUB auctions. Data regard all auctions organized during 2007, 2008, 
2009 and part of 2010 by these web sites. We collected detailed information concerning the auctions: the value of the goods, the 
cost of the fee, the maximum bid amount, the duration of the auction or eventually the required number of bids. We report in 
Table [ST] some of these quantities calculated for our data sets. We were able also to keep track, for all data sets, of each single 
bid, getting information about its value, the time when it was made and the agent who made it. Data sets were anonymised and 
can be downloaded at the page f ilrad . homelinux . org. In the following, we focus our analysis mainly on the data sets 
UBH and BM since, given their size, allow to perform much better statistics. 
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C. Analysis 
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Figure S2: UBH data set. (a) In the top panel, the time series of bid values performed by agent u = 1632 in auction a = 19 is shown. In 
particular we zoom into a region where "abnormal" movements are present. In the bottom panel, we plot the time series of the length of her 
jumps dt = \bt+i — bt\. (b) Same plots as those appearing in panel a, but for the cleaned version of the time series. The indices a and u refer 
to the anonymised version of the UBH data set. 



Fixed an auction a and an agent u, we consider the temporal series of her bids, whose total number is denoted by Ta^w This 
series is basically a list of integers 6i, 62, . . . , 6t, where we have suppressed the indices a and u for shortness of notation. Their 
value is defined over the interval [1, M]. All bs are different each other since no agent bids on a certain value more than once. 
A typical example of these time series is reported in the top panel of Figure [S2k. We calculate the gap or difference between 
subsequent bid values and indicate it with dt = |6t+i — bt\. Given a list of T bid values we can in fact extract a list of T — 1 
differences between consecutive bids. In the bottom panel of Figure S2 i, we plot the time series dt for the same agent whose 
bid time series is plotted in the upper panel. 
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D. Cleaning the data sets 

Generally, more "professional" agents perform random searches followed by systematic coverages of intervals. A typical ex- 



ample is shown in the inset of the upper panel of Figure [S2fi . Here a zoom of the time series appearing in the main plot is 
reported. Systematic coverages are performed by selecting a range of bid values and then placing a bid on each single value in 
that interval. This is an opportunity offered by the web site hosting auctions. Each bid in this case is characterized by the same 
time stamp. Such occurrence is likely for agents who make a significant number of bids, but becomes less relevant for agents 
who invest relatively small amount of money. We cleaned data by removing all pieces of the time series corresponding to this 
"abnormal" behavior. It should be remarked that the gaps between consecutive bids which can be measured in these regions are 
in the majority of the cases equal to one. Including these regions will influence only gaps equal to unity by overestimating their 
presence. We decided to remove such systematic coverages in order to focus our attention only on "normal" bidding strategies. 
The result, after the cleaning procedure of time series shown in Figure |S2k, is reported in the upper panel of Figure |S2b. The 



gaps between consecutive bids in the cleaned time series are reported in the bottom panel of Figure [S2| 



As additional information, in Figure S3 we measure the number of agents N (p) performing a ratio p of bids made using a 





0.2 0.4 0.6 0.8 



0.2 0.4 0.6 0.8 1 



Figure S3: Number of agents N (p) performing "normal" search strategies at rate p. 



normal strategy (i.e., after removing systematic coverages) and the total number of bids placed. 

In the following analysis, we consider only cleaned time series, where systematic coverages have been deleted. 
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E. Statistics of the length of the jumps 
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Figure S4: (a) Probability distribution function P (d) calculated over all agents and auctions in the data sets UBH (orange circles), LB (gray 
squares) and BM (turquoise diamonds). Dashed lines stand for best power-law fits (least square). We find a = 1.54(2) [black], a = 1.54(3) 
[red] and a = 1.54(5) [blue], (b) Same as in panel a, but for cleaned time series. Dashed lines have slopes a = 1.49(2) [black], a = 1.52(2) 
[red] and a = 1.51(5) [blue]. The data of this figure also reported in Figure 3A of the main text. Curves calculated for LB and BM data sets 
have been vertically shifted for clarity. 



We measure the probability distribution function (pdf) P (d) of the difference between subsequent bids made by single agents 
in single auctions. Global (i.e., aggregated over all agents and auctions) pdfs of both data sets are plotted in Figure [S?) For 
agents with a sufficient number of bids in the same auction, we also compute individual pdfs. Some examples are reported in 



the various panels of Figures [S5][S6] and |S7]for UBH data set and in Figures [S8]|S9] and |SlQ| for BIM data set. In each panel of 
these figures, we explicitly indicate the id of the auction a and the id of the agent u as they appear in our anonymised version of 
the data sets. In all cases, we find a behavior compatible with 



P {d) - . 



(SI) 



The search strategy adopted by agents is therefore given by Levy flights with characteristic exponent a. The best fits with 
power-laws are plotted in Figures [SS] [56] [ST] [SS] [S9| and |SlQ| with black dashed lines and the value of a, plus the associated 
error, corresponding to the best fit is reported at the top of each panel. 

In order to calculate a pdf, we divide the range of possible values of d in bins equally spaced on the logarithmic scale. We then 
drawn the pdf by associating to each bin the number of ds falling in that bin divided by the number of integers that could enter in 
the bin (this because d can assume only integer values). Everything is then normalized by simply dividing by the total number of 
points. We compute the best power-law exponent by performing a linear least square fit in double logarithmic scale. Despite the 
binning procedure may introduce a certain amount of arbitrariness in the evaluation of the pdfs, we checked the consistency of 
our results by varying the number of bins. IMoreover, we additionally make use of a different fit method (maximum likelihood) 
about which we will discuss later. 
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a=l u = 23 a =1.6(1) 



a = 5 u=128 a =1.6(1) 





a =19 u=l a =1.3(1) 



a =19 u=134 a =1.2(1) 
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Figure S5: UBH data set. Probability distribution function P (d) measured for agent u in auction a. We show several P {d)s for different pairs 
u and a. Dashed lines denote the best power-law fit (least square) P (d) ^ d~^ obtained for the data. 
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19 u=1601 a =1.4(1) 



1=19 u= 1632 a =1.4(1) 
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a =19 u= 1642 a =1.3(1) 
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Figure S6: UBH data set. Same as Figure [S5] 



:73 u= 1640 a =1.5(2) 



1 = 91 u=1715 a =1.5(2) 





d d 

Figure S7: UBH data set. Same as Figure[S5]and |S6] 
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a =105 u = 28 a =1.5(1) 



a =122 u=15 a =1.5(1) 





a =125 u = 28 a =1.4(1) 



a = 293 u = 3503 a =1.5(1) 





a = 300 u = 550 a =1.5(1) 



a = 327 u=150 a =1.6(1) 



CLh 10 





Figure S8: BM data set. Probability distribution function P (d) measured for agent u in auction a. We show several P {d)s for different pairs 
u and a. 
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a = 50 u = 48 a =1.5(1) 



a = 55 u=15 a =1.3(1) 





a = 61 u = 213 a =1.3(1) 



a = 61 u = 36 a =1.3(1) 





a = 62 u=136 a =1.6(1) 



a = 69 u = 28 a =1.5(1) 





Figure S9: BM data set. Same as Figure [S8] 
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a = 89 u = 48 a =1.5(1) 



a = 89 u = 72 a =1.4(1) 





a = 91 u=15 a =1.4(1) 



a = 94 u=15 a =1.6(1) 





Figure SIO: BM data set. Same as Figure [S8] and |S9] 



19 



F. Independence of the direction of tlie jumps 
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Figure Sll: (a) UBH data set. Probability distribution function P (d) calculated over all agents and auctions. We have separated positive from 
negative variations. The measured decay exponents of the best fits with power-laws (dashed lines) are a = 1.53(2) [black] and a = 1.54(2) 
[red], respectively. The curve corresponding to negative variation has been vertically shifted for clarity. This figure appears also in Fig. 3B of 
the main text, (b) BM data set. here the best power-law fit are a — 1.61(8) [black] and a — 1.45(4) [red]. 



Since the rules of the auction naturally bring agents to move towards low bid values, it is important to stress any eventual 
difference between the statistics associated with the length of gaps between consecutive bid values. We aggregated data from all 
agents and auctions and separate positive (i.e., > ht) from negative (i.e., < ht) variations. In Figure Sll the pdfs of 
positive and negative variations are plotted together. As one may notice, there is not a significant difference between them and 
both show a clear power-law decay with compatible exponents. 
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G. Independence of the agents' activity 
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Figure S12: (a) UBH data set. Probability distribution function P (T) of the number of bids T performed by single agents in single auctions. 
We calculate P (T) on the original data set (orange circles) and on the cleaned version of the same data set (gray squares). In both cases, the 
distribution scales power-like with a decay exponent equal to 2.2(2). (b) BM data set. The exponent of the best fit is 2.0(4) [dashed line]. 



The evidence of Levy flights for single agents in single auctions can be directly verified only for agents with a sufficient number 
of bids T in the same auction. For values of T smaller than 50 is practically impossible to construct the histogram P (d) and 
therefore no exponent can be measured. Unfortunately, this situation is very frequent in our data sets. We measure the number 
of bids made each agent in each auction and plot the pdf of the number of bids in a single auction in Figure |S12[ The level of 
activity is quite heterogeneous and decays power-like with an exponent close to 2.2. For completeness, we measure the same 
pdf for both original and cleaned data sets without noticing appreciable differences. 
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Figure S13: (a) UBH data set. Probability distribution function P (d) of the length d of the jumps performed by agents in the bid space. All 
auctions have been aggregated together. Different curves correspond to agents with different levels of activity. Their activity is measured as 
the number of bids made in the same auction. We divide the population into four subsets: T < 10 (orange circles), 10 < T < 40 (gray 
squares), 40 < T < 200 (blue diamonds) and T > 200 (violet triangles). Dashed lines represent the best power-law fits. The value of the 
measured exponents are: a = 1.55(2) (black), a = 1.62(4) (red), a = 1.53(2) (blue) and a = 1.54(3) (violet). Curves have been vertically 
shifted for clarity. This figure is also reported in Fig. 3C of the main text, (b) BM data set. Best power-law fits (dashed lines) have exponents: 
a = 1.43(4) (black), a = 1.44(5) (red), a = 1.7(1) (blue) and a = 1.7(1) (violet). 



We divide the population in different ranges of activity. We aggregate the length of the jumps performed by all agents in a given 
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bin and measure the resulting P (d). The results of this analysis are reported in Figure S13 Independently of the activity level, 
the aggregated pdfs decay power-like and have similar exponents a. This means that the presence of Levy flights is typical for 
every agent independently of how many bids the agent has made. 
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H. Independence of the bidding time 
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Figure S14: (a) UBH data set. Probability distribution function P (d) of the length d of the flights performed by agents in the bid space. All 
auctions have been aggregated together. Different curves correspond to different periods of activity for agents: orange circles correspond to 
the initial guess made by agents; bids gaps with t < 2 (gray squares) and t < 5 (blue diamonds) aggregate the data corresponding to the early 
activity of agents; t > T — 2 (violet up triangles) and t — T — 1 (green down triangles) corresponds to the jumps made by agents at the 
end of their own activity. Dashed lines have been obtained as best power-law fits with data points. The value of the measured exponents are: 
a = 1.55(4) (black), a = 1.54(3) (red), a = 1.59(4) (blue), a = 1.49(5) (violet) and a = 1.54(6) (green). Curves have been vertically 
shifted for clarity. This figure is also reported in Fig. 2d of the main text, (b) BM data set. Best power-law fits (dashed lines) have exponents: 
a = 1.3(1) (black), a = 1.50(5) (red), a = 1.47(4) (blue), a = 1.43(6) (violet) and a = 1.35(6) (green). For the initial bid 6i, we compute 
the length of the jump as (i = M — 6i + 1, with M being the maximal bid value in the HUB auctions. 



Another fundamental point is to understand whether the Levy flight strategy is emergent or a priori given. We test these 
hypotheses by measuring the pdfs P {d) corresponding to a certain range during the activity of the agents. The results are 



reported in Figure S14 We consider ranges of activity periods corresponding to t < 2, t < 5,t > T — 2 and t = T — 1. 
Additionally, we consider the distribution of the initial bid values (i.e., the first bid made by all agents in all auctions). In every 
case, we are able to fit the curves with power-laws and the resulting exponents are compatible each other. We can effectively 
conclude that the strategy to adopt a Levy flight is not an emergent property induced by the evolution of the auction. Instead the 
strategy to follow Levy flights is intrinsically present, in each agent, during the whole duration of the auction. 
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I. Independence between jumps 




T = 5 T=10 




Figure S15: UBH data set. Number of agents N (r) whose bid gaps at position t and t + r have Pearson's correlation coefficient equal to r. 
Black curves are measured on real data, while the red ones are calculated over a reshuffled version of the same data. The reshuffling is made 
by randomly exchange pairs of entries in the time series of the time gaps with the only prescription that the sum of them is not lower than one 
and not larger than M. We consider different values of r. In each plot only agents with at least 10 + r bids in the same auction are considered. 



We further study the correlations between jumps. Given an agent and an auction, we consider the list of all her jumps 
(ii, (i2, . . . , <^B-i and calculate the Pearson's correlation coefficient 



{{dt - lit) {dt+r - l^t+r)) 



(S2) 



where (•) stands for the average over the entire time series (i.e., over all values of t from 1 to T — 1 — r). /i^ = {dt) and jit+r 



{dt-\-r) are the average values of the bid gaps along the time series, while at = \/ {d^) — {dt) and cjt+T = \J (<^t+r) ~ {^t+r) 
are the respective standard deviations. We measure such coefficient for every agent who has performed at least 10 + r bids 



in the same auction and show the number of agents N (r) with given value of r in Figures S15 and S16 The same quantity 
is also calculated for a randomized version of the time series, where bid gaps are randomly reshuffled with the only constraint 
that their partial sum cannot never be smaller than one and larger than M. We consider several values of r. As one can clearly 
notice, subsequent gaps (i.e., r = 1) are slightly correlated. Such correlation, becomes negligible when r grows and already for 
r = 2, A/" (r) is negligible. For r = 5 and r = 10, the curves corresponding to the original time series and those obtained over 
randomly reshuffled time series are almost identical. 

Such results show that agents perform almost uncorrelated Levy flights. Once an agent makes a jump, the length of this jump 
is slightly correlated with the one of the jump made before. However, after few jumps there is not longer memory of what 
happened before. In good approximation, the walk of the agent in the bid space can be therefore modeled as the one followed 




random walker performing uncorrelated Levy flights. 
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J. Testing the model 



1. Maximum likelihood fit and Goodness of fit 



In this section, we compute the level of significance of our model for the description of real time series. Suppose that a time 
series of T bid values 61 , 62, • • • , describes a realization of our model. Fixed the exponent a, the bound M of the lattice and 
the position ht-i at stage t — 1, the probability that the random walker jumps at ht at stage t is given by the transition matrix of 
Eq. (6) of the main text. The probability or likelihood that the whole sequence was extracted from our model is 



{Qa)b 



The value of a that maximizes the former equation represents the most likely exponent of our model that could have generated 
our particular sequence. In order to find its maximum, it is convenient to take the logarithm of both sides and write the log- 
likelihood 



T 

£(61,62,. . .,bT\a) = Y^ In \ {Qa)b,_,,bt = X^l^l^t - ^t-i \ - X^ln [f^h-i i^)] , 



^=1 



(S3) 



where we set bo = 0. The value a' at which the maximum of Eq. (S3) occurs can be estimated numerically. 
a' is the best exponent fitting the data in the hypothesis that they were produced according to our model. The significance level 
of the model for the description of the data can be calculated by estimating the p- value associated with our measurement. In this 
respect, we first compute the distance between the theoretical distribution of the jump lengths 



Pa' (d) 



d i,j 



and the one obtained from our data P (d) by calculating 



rjdata = max|Pc,/ {> dt) 



P{>dt 



where P {> d) = Xlg>d P (q) the cumulative distribution of the jump lengths. Notice that the distance between cumulative 
distributions is the same as the one adopted in the Kolmogorov-Smirnov test. We then generate artificial time series of length T 
from our model with exponent and compute their distance r] with respect to the theoretical distribution. The p- value is finally 
determined by the relative number of times in which we observe r] > r]data- 

Synthetic time series generated according to our model can be additionally used for the determination of the error associated 
to the estimation of a'. The error associated to a' is the standard deviation of the best exponents estimated, with maximum 
likelihood, for the synthetic data sets. 

A graphical comparison between the exponents a (least square method) and a' (maximum likelihood method) is presented in 
Fig ure |S17| In general, the two methods produce consistent results. For completeness, we list the results obtained in Tables [S2 
and[S3l 

The p- values show also a general goodness of our model for the description of the data. Sometimes however, the value of p 



is very small. This could be explained in a simple manner. In Figure [ST8] we plot for example the statistical test performed over 

13 000. However, this value 
2 000 as the sudden 



the same time series appearing in Figure S2 In that auction, the maximum bid amount was M 
of M does not correspond to the "effective" bound felt by the agent. This bound seems to be around M 
drop of P (> d) would suggest. By setting M = 2 200 and running again the statistical test, as we did in Figure S18 3, we see 
clearly that the curve predicted by our model and the one measured on real data are in very good agreement. 
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Figure S17: We consider only agents who have performed at least T = 50 bids in a single auction in UBH data set and those with at least 
T = 100 in the BM data set. Our data set offers 39 agents that satisfy this constraint in the UBH data set and 52 in the BM data set. Panels 
a (UBH) and b (BM) show the scatter plot ol versus a for the best power-law exponents estimated by using maximum likelihood and least 
square methods, respectively. The agreement between the two measurements is good as demonstrated by the fact that the majority of the points 
fall on the diagonal (red line). 
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Table S2: UBH data set. Each row corresponds to one of the 39 agents who have bid at least 50 times in the same auction. We report the id 
of the auction a, the id of the agent u, the exponent a calculated with the least square method, the exponent a calculated with the maximum 
likelihood method and the p- value. Entries with low p- values are marked with * . In the 77% of the cases we find a p-value larger than 0, which 
indicates that our model well describe the time series. 
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Table S3: BM data set. Each row corresponds to one of the 52 agents who have bid at least 100 times in the same auction. We report the id 
of the auction a, the id of the agent u, the exponent a calculated with the least square method, the exponent a' calculated with the maximum 
likelihood method and the p- value. Entries with low p- values are marked with * . In the 79% of the cases we find a p- value larger than 0, which 
indicates that our model well describe the time series. 
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Figure SI 8: Example of the maximum likelihood method for the determination of the exponent a' of the Levy flight. We consider the time 
series of agent 1632 in auction 19 of the UBH data set (the same appearing in Figure [S2|. In panel a, we set M = 13 000 which is the value 
of the maximum bid amount that was allowed in the auction. The best exponent a — 1.2(1) is obtained by looking at the maximum of the 
log-likelihood as it is shown in the inset. The comparison between the cumulative distribution of the jump lengths expected from the model 
(red line) and the one calculated over the time series (black) do not well agree. It seems that the "effective" bound is smaller than the real 
one. In panel b, we set M = 2 200 and consider only bid values smaller than this bound. On the new time series, we perform a maximum 
likelihood fit finding — 1.1(1). The theoretical expectation (red line) and the one obtained from the time series (black line) are now very 
similar yielding a p- value equal to 0.1. 
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2. Another maximum likelihood fit 



Since the upper-bound M is agent dependent, we perform an additional analysis where the upper bound M is not directly taken 
from the data, but used as a parameter for the fit. For each agent, we let the parameter M vary only in the range for which at 
least the 90% of the bids values are below M. Indicate with T the number of bids below the threshold M. We then find a' 
identifying the maximum of the likelihood function and cal culate the p- value as described so far. We consider the best value of 



M as the one which maximizes the product T x p. Figures 



S19 



S20 



and 



Figures [S5| [S6l and [ST] for UBH data set, while Figure s |S22[|S23| and [S24[are the analogous of Figures [S8||S9] and [STO| for BM 



S21 



report the best fit for the same agents analyzed in 



data set. For completeness, we report in Tables [S4] and |S5] the results obtained with the maximum likelihood fit where M is used 
as parameter of the fit. The best exponents are still consistent with those reported in Tables S2 and S3 but the p- values result 
much increased. 
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Table 84: UBH data set. Each row corresponds to one of the 39 agents who have bid at least 50 times in the same auction. We report the id 
of the auction a, the id of the agent u, the exponent a calculated with the least square method, the exponent a calculated with the maximum 
likelihood method, the best value of the upper bound M and the p- value. Entries with low p- values are marked with * . In the 90% of the cases 
we find a p- value larger than 0, which indicates that our model well describe the time series. 




Figure S19: UBH data set. Cumulative distribution function P {> d) measured for agent u in auction a (red full line) compared with the 
theoretical distribution (black dashed line). We show several P (> c/)s for different pairs u and a. We report also the best value of the upper 
bound M and the p- value associated with our fit. 
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a = 73 u= 1640 a =1.4(1) M= 1920 p-val = 0.27 



a = 91 u=1715 a =1.3(1) M = 70 p-val = 0.67 





d d 

Figure S21: UBH data set. Same as Figure [819] and |S20| 
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0.17 


91 


15 


1.4 


;i) 


1.4( 


;i) 


180 


0.87 


92 


36 


1.8 


[2) 


1.3( 


1) 


130 


0.04 


92 


28 


1.3 


;i) 


1.4( 


1) 


910 
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15 
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28 
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98 
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Table S5: BM data set. Each row corresponds to one of the 52 agents who have bid at least 100 times in the same auction. We report the id 
of the auction a, the id of the agent u, the exponent a calculated with the least square method, the exponent a calculated with the maximum 
likelihood method, the best value of the upper bound M and the p- value. Entries with low p- values are marked with * . In the 92% of the cases 
we find a p- value larger than 0, which indicates that our model well describe the time series. 



33 



a =105 u = 28 a =1.5(1) M=1310 p-val = 0.59 a =122 u=15 a =1.4(1) M = 470 p-val = 0.95 




d d 

Figure S22: BM data set. Cumulative distribution function P {> d) measured for agent u in auction a (red full line) compared with the 
theoretical distribution (black dashed line). We show several P (> c/)s for different pairs u and a. We report also the best value of the upper 
bound M and the p- value associated with our fit. 



a = 89 u = 48 a =1.3(1) M = 230 p-val = 0.42 a = 89 u = 72 a =1.3(1) M = 330 p-val = 0.30 




Figure S24: BM data set. Same as Figure |S22] and |S23| 
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3. Probability distribution of the Levy flight exponents 




a a 



Figure S25: The pdf g (a) of the exponents calculated with the three fitting methods. The pdfs corresponding to the maximum likelihood fit, 
where M is a parameter of the fit, are the same as those appearing in Figs. 2D and 2E of the main text. The distributions are characterized 
by the following values of the mode at, average {a) and variance a. UBH data set: for least square fit we have at = 1.34, (a) = 1.40, 
a = 0.26; for maximum likelihood fit we have at = 1.28, (a) — 1.29, a — 0.20; for maximum likelihood fit with additional fitting parameter 
M we have ah = 1.21, {a) = 1.26, a = 0.23. BM data set: for least square fit we have ah = 1.46, {a) = 1.42, a = 0.27; for maximum 
likelihood fit we have ah — 1.34, {a) — 1.37, a — 0.21; for maximum likelihood fit with additional fitting parameter M we have ah — 1.35, 
(a) = 1.36, cr = 0.23. 



As a final result, we compute the distribution g {a) of the exponents measured for single agents in single auctions. We still 
consider only agents who have performed at least T = 50 bids in a single auction in UBH data set and those with at least 
T = 100 in the BM data set. Assuming that the best estimation of the exponent of agent u is au and the associated error of the 
measurement is Aa^, the empirical distribution of the exponents is calculated as 

g {a) = ^ e-(«.-)^ / [2 (Aa^f] ^ (S4) 

where C = J da g- / [2(Aq!^) ] ^j^^ proper normalization constant. The resulting pdfs are reported in 

Fig.[S25| 



